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The  work  reported  here  is  taken  verbatim  from  a proposal  for  the  study  and  de- 
sign of  a computer  aided  sonar  echo  classification  system.  In  the  proposal,  the 
discussion  of  signal  processing  and  waveform  design  was  felt  to  be  of  sufficient 
general  nature  to  be  useful  to  others  working  in  related  fields.  Sid  Applebaum, 
who  has  since  terminated  his  services  with  the  Company,  is  the  sole  author  of  this 
work. 


TIS  R64EMH11 

SIGNAL  PROCESSING  AND  WAVEFORM  DESIGN 


The  intent  here  is  to  discuss  some  theoretical  aspects  of  the  signal  process- 
ing and  waveform  design  areas.  In  addition,  some  recent  findings  on  "noise"  wave- 
forms shall  be  presented.  Noise  or  random  waveforms  are  especially  attractive  for 
target  classification  because  they  combine  the  attributes  of  long  duration  and  wide 
bandwidth,  thus  offering  the  possibility  of  simultaneous  resolution  in  both  range 
and  doppler.  Noise  waveforms  are  attractive,  too,  because  techniques  have  been  de- 
veloped for  efficiently  generating  and  processing  such  waveforms. 

It  is  convenient  for  this  discussion  to  group  reflectors  of  sound  energy  into 
two  classes  - point  targets  and  extended  targets.  A point  target  is  one  whose 
physical  dimensions  are  small  compared  to  a wavelength  of  the  transmitted  carrier 
frequency,  while  an  extended  target  has  large  dimensions  compared  to  a wavelength. 

If  acoustic  radiation  is  incident  on  an  extended  target  having  a complicated  shape, 
it  is  fcund  that  the  main  contribution  to  the  scattered  field  is  due  to  "specular" 
points,  or  patches,  occupying  a relatively  small  part  of  the  Illuminated  surface  of 
the  object.  An  analogous  phenomenon  of  comon  experience  occurs  in  the  "spotty" 
reflection  of  the  sun  or  moon  from  a wavy  water  surface. 

The  number  and  location  of  the  specular  points  observed  depends  upon  the  tar- 
get's shape  and  size,  and  on  its  orientation  with  respect  to  the  observer.  Objects 
of  simple  geometric  shape  will  exhibit  a small  number  of  specular  points.  A sphere, 
for  example,  has  only  one  specular  point  and,  therefore,  will  look  like  a point  tar- 
get even  though  it  may  be  large  compared  to  a wavelength. 

Extended  targets  may,  therefore,  be  considered  to  be  a group  or  cluster  of 
point  targets.  The  distribution  of  these  point  targets  in  range  and  amplitude  will 
be  characteristic  of  the  target  shape,  and,  hence,  may  be  used  to  recognize  or  clas- 
sify the  target.  The  function  of  signal  processing  in  classification  is  to  obtain 
this  distribution  accurately  and  in  as  much  detail  as  the  transmitted  waveform  per- 
mits. The  waveform  should  be  capable  of  providing  high  resolution  information, 
preferably  enough  to  individually  resolve  at  least  a few  of  the  target's  major 
specular  points.  If  a sonar  system  has  poor  range  resolution,  all  of  the  target's 
specular  points  may  appear  as  a single  point  target,  and  there  will  be  little  infor- 
mation on  which  to  base  a classification  decision. 
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Target  motion  can  be  of  considerable  help  in  classification.  If  a rigid  ex- 
tended target  is  moving,  the  signals  from  all  of  its  specular  points  will  display 
the  same  doppler  shift.  Doppler  shift  can  thus  serve  as  a tag  to  identify  the 
specular  point  returns  from  a single  extended  target.  With  sufficient  doppler 
resolution,  the  group  of  specular  points  belonging  to  one  moving  extended  target 
can  be  separated  from  the  group  belonging  to  another  stationary  target,  even  though 
both  targets  may  be  at  the  same  range.  Doppler  information  can  also  help  recognize 
"nonrlgld"  targets  like  a school  of  fish,  on  the  basis  of  the  amount  of  doppler 
spread  of  the  returned  echoes. 

Since  all  targets  can  be  considered  as  groups  of  point  targets,  the  waveform 
design  and  signal  processing  portions  of  the  classification  problem  reduces,  in 
effect,  to  the  problem  of  detecting,  resolving,  and  estimating  the  parameters  of 
point  targets.  A great  deal  of  work  has  been  done  in  this  area.  A brief  review  of 
the  major  theoretical  results,  and  their  application  to  the  problem  of  sonar  echo 
classification  is  presented  in  the  remainder  of  this  section. 

A.  OPTIMUM  SIGNAL  PROCESSING 

A discussion  of  signal  processing  and  waveforms  is  necessarily  mathematical. 
To  save  time  and  space  complex  signal  notation  shall  be  used.*  In  complex  notation, 
the  transmitted  signal  may  be  represented  as 

(2  Et)1/2  u(t)  exp  j 2n  fc  t (1) 

where  u(t),  a complex  waveform  is  normalized  so  that 

/ | u(t) | 2 dt  - 1 (2) 

E^  - transmitted  energy,  and 
f£  * the  carrier  frequency. 

(Unless  otherwise  specified,  the  limits  on  all  integrals  are  to  be  taken  from  -o  to 

If  there  is  a stationary  point  target  present  at  a round  trip  range  delay 
of  r seconds,  the  signal  received  from  the  target  may  be  represented  by 

s(t)  exp  J 2n  f£  (t  - r)  (3) 

*See  P.M.  Woodward,  "Probability  and  Information  Theory  with  Applications  to  Radar", 
Pergamon,  London  1953. 
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where 


s(t)  - (2  Eg)1/2  u (t  - r)  (4) 

Eg  is  the  received  signal  energy.  Eg  depends  upon  the  transmitted  energy, 
on  the  cross  section  of  the  target,  and  on  attenuation  factors  that  are  functions 
of  range.  It  has  been  assumed  that  the  signal  bandwidth  is  sufficiently  narrow  so 
that  dispersion  in  the  medium  produces  negligible  distortion  of  the  signal  waveform, 
u(t)  . 

If  the  point  target  is  moving  with  a constant,  positive  range  rate,  the  re- 
ceived signal  will  exhibit  a corresponding  doppler  frequency  shift,  say  -f^.  The 
complex  envelope  of  the  received  signal  will  then  be  of  the  form 

s(t)  - (2  E,,)1^2  u(t  - r)  exp  -j  2 n (t  - r)  (5) 

It  is  apparent  that  the  signal  from  any  point  target  is  completely  speci- 
fied by  three  parameters  - energy,  range  delay,  and  doppler  frequency  shift.  All 
of  the  classification  information  available  from  a resolved  point  target  echo  is 
contained  in  the  parameters  Eg,  r,  and  f<j,  or  functions  of  these  three  parameters. 

The  information  available  from  a cluster  of  unresolved  signal  echoes  is  completely 
defined  by  the  distribution  of  received  energy  in  range  delay  and  doppler  frequency. 
Both  types  of  information  can  be  provided  by  a signal  processor  that  estimates  the 
received  energy  Eg  continuously  as  a function  of  both  r and  f^. 

The  type  of  signal  processing  and  waveform  to  use  for  optimum  parameter 
estimation  depends  on  the  environment  in  which  the  system  operates.  In  addition  to 
the  signal  from  the  target  of  interest,  the  receiver  will  pick  up  signals  from  other 
sources,  such  as  neighboring  targets,  reverberation,  multipath  signals,  and,  of 
course,  random  thermal -type  noise.  These  environmental  factors  determine  the  threshold 
detectability  of  a target  and  the  accuracy  with  which  its  parameters  can  be  estimated. 

The  situation  in  which  random  noise  is  the  major  environmental  factor  shall 
be  considered  first.  The  signal  processing  part  of  the  classification  problem  then 
reduces  to  the  well-known  problem  of  estimating  the  parameters  of  a signal  in  addi- 
tive Gaussian  noise.  The  waveform  received  in  this  case  contains  two  components, 
and,  omitting  the  carrier  frequency  factors,  may  be  represented  as 

v (t ) - s ( t ) + n(t)  (6) 

where  n(t)  “ complex  noise  modulation. 
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It  has  been  shown  in  the  literature  that  the  optimum  signal  processor  for 
this  situation  is  one  that  cross-correlates  the  received  waveform,  v(t),  with  a 
stored  replica  of  u(t),  the  transmitted  waveform.  The  cross-correlation  must  be 
performed  for  all  range  and  doppler  shifts  of  interest.  The  optimum  signal  proc- 
essor can  be  defined  more  precisely  as  one  that  generates  a two-dimensional  array 
of  data,  x(r,  ?j)  where 


*(r,  fd) 


^ j 2n  f . t 

/ v(t)  u*  (t  - r)  e dt 


'?) 


The  cross-correlation  can  be  performed  "actively"  by  storing  u(t)  at  the 
receiver,  multiplying  by  u(t  - r)  (with  r varying  over  the  range  of  interest),  and 
then  integrating  the  product  in  some  form  of  spectrum  analyzer.  The  Deltic-VICI 
technique  is  one  method  for  accomplishing  this.  The  cross-correlation  can  also  be 
performed  "passively"  in  a bank  of  "matched  filters".  In  this  case  the  receiver 
consists  of  an  array  of  filters  with  each  one  conjugate-matched  to  a frequency- 
shifted  equivalent  of  the  transmitted  waveform,  u(t).  The  transfer  functior  of 
each  filter  will  be  of  the  form 

u.  (f  - f.  - fd> 


where  U(f)  is  the  Fourier  transform  of  u(t),  that  is 

U(f)  - / u(t)  e'j  2n  £ t dt  (8) 

The  number  of  filters  required  in  a "matched-filter"  receiver  is  approximately  equal 
to  the  total  doppler  surveillance  band  divided  by  the  system  doppler  resolution, 
which  depends  on  the  duration  of  the  waveform,  u(t) . In  the  matched-filter  receiver 
the  doppler  shift  of  the  signal  is  determined  by  noting  which  filter  has  the  largest 
output,  range  is  determined  by  the  time  at  which  the  peak  output  occurs,  and  energy 
is  determined  by  measuring  the  amplitude  of  the  peak  output. 

The  cross-correlating  or  matched-filter  receiver  is  characterized  by  the 
property  that  it  maximizes  the  ratio  of  peak  signal  output  to  rms  noise  output.  The 
mean  square  noise  output  will  be  the  same  for  all  values  of  r,  fd  and  is  given  by 

Xn  " No  / lU<f)!2  df  <9> 

- Nq  / | u(t)| 2 dt 

where  NQ  is  the  (one-sided)  power  density  of  the  noise  (watts  per  cycle). 
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The  "signal"  portion  of  the  signal  processor  output  will  be 


j 2 tt  f t 

Xg(r,  ?d)  * / s(t)  u*(t  - r)  e dt 


(10) 


If  s(t)  is  the  signal  returned  by  a point  target  at  a range  delay  r,  moving 
so  as  to  produce  a dcppler  shift  of  -fd,  Equation  (10)  becomes 


:(r, 


(2  Eg) 1/2  / u (t 


r)  u*(t  - r)  e 


-J  2,  (fd-£d) 


dt  (11) 


or 


XsCr, 


fd>  ’ (2  ES)1/2  Au  <r  * 


fd> 


(12) 


where 


A (t,  v)  - / u(t)  u*(t  + r)  e"J  2n  Vt  dt  (13) 

u 7 

A ^ 

The  peak  signal  output  will  occur  when  r ■ r,  and  fd  ■ fd)  at  which  point, 
Xs(r,  fd)  - (2  Es)1/2  *u<0,0)  (14) 

-(2  «,)l/2 


since  Ay (0,0)  * 1,  from  Equations  (2)  and  (13). 


The  ratio  of  peak  signal  output  to  rms  noise  output  from  Equations  (9)  and 
(14)  is  therefore 


Peak  Signal 
rms  Noise 


(2 


V 


1/2 


1/2 


(15) 


This  is  an  important  fundamental  result  showing  that  the  output  signal -to- 
noise  ratio  depends  only  on  the  energy  in  the  received  waveform  and  not  at  all  on 
its  shape.  Since  the  peak  power  that  can  be  transmitted  is  limited  by  cavitation, 
high  slgnal-to-ncise  ratios  can  only  be  obtained  by  using  waveforms  of  long  duration. 
Long  waveforms,  however,  have  poor  range  resolution  unless  some  form  of  coding  or 
modulation  is  used.  For  target  classification,  it  is  desirable  to  have  both  high 
signal -to-noise  ratio  and  good  range  resolution.  Long  coded  waveforms  must  there- 
fore be  considered,  that  is,  waveforms  that  have  a large  bandwidth-time  product. 
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The  resolution  capabilities  of  a waveform  in  a matched  filter  system  are  determined 
by  the  function  Ay^v)  defined  in  Equation  (13). 

The  function  Ay  is  called  the  ambiguity  function  of  the  waveform  u(t);  it 
gives  the  distribution  of  "interference"  or  "crcss-talk"  due  to  a single  point  tar- 
get. The  ambiguity  function  is  a sort  of  aperture  or  window  with  which  targets 
distributed  in  range  and  velocity  are  viewed.  From  the  definition,  Equation  (13), 
the  following  important  relations  may  be  derived 


au(t,v)  - / U(f  + v)  U*(f)  e‘J  2n  fT  df 

(16) 

|Au(t,v)|2  < |Au(0,0)|2 

(17) 

Iau(-t,-v)|2  - |Au(t, v) | 2 

(18) 

If  I Au (t , v) | 2 dx  dv  - | A(0,0) 1 2 - 1 

(19) 

Equation  (16)  expresses  the  ambiguity  function  in  terms  of  U(f),  the 
Fourier  transform  of  u(t) . Equation  (17)  shows  that  the  ambiguity  function  attains 
its  peak  value  at  the  "target  position"  (t,v  ■ 0,0),  verifying,  in  part,  the  opti- 
mality of  matched  filtering.  Equation  (18)  reveals  an  Important  reciprocal  property 
of  the  ambiguity  function;  it  can  be  shown  from  Equation  (18)  that  the  cross-talk 
power  at  the  point  t,v,  due  to  a target  at  0,0,  equals  the  cross-talk  power  at  0,0, 
due  to  an  equal  size  target  at  the  point  t,v.  Equation  (19)  is  known  as  the  Uncer- 
tainty Principle;  it  says  that  the  total  integrated  cross-talk  power  generated  by  a 
point  target  is  a constant  independent  of  the  shape  of  the  waveform  u(t). 

Although  the  total  cross-talk  power  or  ambiguity  is  invariant  with  the 
choice  of  waveform,  u(t),  the  manner  in  which  it  is  distributed  in  range  and  doppler 
depends  very  much  on  the  waveform  coding.  Ideally  it  would  be  preferred  for  A (t,  v) 
to  consist  of  a single  sharp  central  spike  with  the  unavoidable  excess  ambiguity 
uniformly  distributed  as  low-level  "side  lobes".  Such  an  ambiguity  function  has 
been  referred  to  as  a "thumb-tack"  function.  The  dimensions  of  the  central  "spike" 
of  an  ambiguity  function  are  approximately  1/B  in  range  (t),  and  1/T  in  doppler  (v), 
where  B and  T are  the  bandwidth  and  time  duration  of  the  waveform,  respectively. 
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The  response  of  a matched  'liter  receiver  to  an  input  signal  pulse  is 
shown  in  Figure  1.  The  array  of  outputs  depicted  defines  the  ambiguity  function  of 
the  pulse  waveform.  The  nature  of  the  ambiguity  functions  of  a number  of  different 
waveforms  is  indicated  in  Figure  2,  In  the  form  of  "ambiguity  diagrams".  In  these 
diagrams,  only  the  regions  of  high  ambiguity  are  shown.  The  side  lobes  of  the 
ambiguity  functions  will  spread  over  2T  in  range  (t)  and  2B  in  doppler  (v) . The 
ambiguity  diagrams  of  a short  and  long  unmodulated  pulse  are  given  in  Figures  2a 
and  2b,  respectively.  The  short  pulse  spreads  its  ambiguity  in  doppler  while  the 
long  pulse  ambiguity  is  spread  in  range.  As  8 result,  the  short  pulse  has  good 
range  resolution  but  poor  doppler  resolution,  whereas  the  converse  is  true  of  the 
long  pulse.  The  ambiguity  diagrams  of  a linear  frequency  modulated  pui&e  are  shown 
in  Figure  2c.  The  region  of  high  ambiguity  in  this  case  has  its  major  axis  along  a 
diagonal  extending  deep  into  the  first  and  third  quadrants.  With  the  linear  FM 
waveform,  resolution  is  poor  for  targets  whose  range-velocity  separations  lie  along 
the  diagonal. 

B.  PERFORMANCE  IN  REVERBERATION 


Up  to  this  point,  random  noise  has  been  considered  as  the  only  environ- 
mental factor  limiting  performance.  Usually,  however,  reverberation  is  at  least  as 
important  as  random  noise  in  determining  sonar  performance.  Reverberation  may  be 
regarded  as  an  extended  target  consisting  of  a large  number  of  small  randomly  dis- 
tributed point  targets  or  scatterers.  Each  point  scatterer  will  produce  a signal 
of  the  form  of  Equation  (12)  at  the  output  of  a matched  filter  receiver.  The  point 
scatterers  are  assumed  to  be  too  densely  distributed  to  be  individually  resolved. 


Because  of  the  random  nature  of  reverberation,  the  average  reverberation 
power  at  the  output  of  a matched  filter  receiver  is  obtained  by  summing  the  indi- 
vidual signals,  power-wise.  If  w(r,f(j)  dr  dfd  is  the  energy  received  from  the  point 
scatterers  located  in  the  range  interval  r,r  + dr  and  producing  doppler  frequency 
shifts  between  fd  and  fd  + dfd,  the  mean  power  out  of  the  matched  filter  receiver 
due  to  reverberation  will  be 


xR(*,*d>  • 2 //  w<r>fd>*!Au<r-r,  VVl2  dr  dfj 


(20) 


In  order  to  proceed  further  with  this  expression,  the  detailed  character  of 
the  function  w(r,f(j)  must  be  known.  For  most  cases  of  Interest  a reasonable  assump- 
tion is  that  w varies  sufficiently  slowly  with  range  so  that  it  may  be  regarded  as  a 
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(A)  SHORT  PULSE  (B)  LONG  PULSE  (C)  LINEAR  FM 


I 

I 

I 

Figure  2.  Ambiguity  Diagrams 
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function  of  fd  only  in  the  integral  shown  in  Equation  (20).  Formally,  this  assumes 

A 

that,  in  the  vicinity  of  r ■ r, 

w(r,fd)  - W(r)  $(fd,r)  (21) 

where  W and  $ are  slowly  varying  functions  of  r,  and  $ is  assumed  to  be  normalized 
so  that 

/ *(fd,r)  dfd  - 1 (22) 


The  function  $ describes  the  distribution  of  the  reverberation  energy  in  the  doppler 
dimension.  Usually  it  will  be  peaked  near  zero  doppler  (assuming  that  any  ship’s 
motion  is  compensated) . Substituting  Equation  (21)  in  Equation  (20) 

X^(r,?d)  " 2 w<r)  If  <Kfd,r)  |Au(r-r,  fd-?d)|2  d*  dfd  (23) 

- 2 W(S)  I ®(fd,r)  Gu(fd-?d)  dfd  (24) 


where 


Gu(v)  « Gu(-v)  - f |A(t,v)|2  dT  (25) 

When  the  distribution  of  reverberation  energy  is  highly  concentrated 
around  zero  doppler,  the  function  $ becomes  a "delta"  function 

® - 5(fd)  (26) 


In  this  case  Equation  (24)  reduces  to 


X2(?,fd)  - 2 W<?)  Gu(1d)  (27) 

These  equations,  (24)  and  (27),  describe  the  distribution  of  reverberation  power  at 
the  output  of  a matched  filter  receiver.  If  there  * a point  target  of  Interest  at 

A A 

r - r,  having  a doppler  shift  fd  ■ fd,  and  returning  a signal  with  energy  Eg,  the 
ratio  of  reverberation  power  to  peak  signal  power  at  the  receiver  output  will  be 


/ «„«<!  - V «„ 


<2» 


for  the  general  case,  and 
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(29) 


W(r) 


W 


when 

• - B(fd) 

Since  both  W and  E§  vary  identically  with  both  the  transmitted  energy  and 
range,  their  ratio  depends  only  upon  the  ratio  of  the  cross  section  density  of  the 
reverberation  scatterers  to  the  cross  section  of  the  target.  The  only  term  in 
Equations  (28)  or  (29)  that  depends  on  the  waveform  u(t)  is  Gy.  This  function  then 
determines  the  ability  of  a matched  filter  system  to  detect  and  classify  targets  in 
a reverberation-limited  environment  and,  therefore,  will  be  called  the  reverberation 

A 

function.  It  varies  with  the  amount  of  target  doppler  shift,  f^.  If  the  reverbera- 
tion is  not  at  zero  doppler,  the  argument  of  the  reverberation  function  Gy  should  be 
interpreted  as  the  difference  between  target  doppler  and  reverberation  doppler.  To 
Indicate  this  more  general  case,  the  notation  Gy(v)  shall  be  used  in  place  of 
(^(f^) . It  can  be  seen  from  Equation  (29)  that  a low  value  of  Gy  indicates  good 
performance  in  reverberation. 

The  reverberation  function  may  be  calculated  from  Equation  (25),  or  from 
one  of  the  following  equivalent  expressions 

Gy(v)  - / |u(f)|2  |U(f  + v)|2  df  (30) 

■ / |au(t,0)|2  e"J  2n  VT  dx  (31) 

If  the  waveform  u(t)  is  not  normalized  according  to  Equation  (2),  then  the 
right-hand  sides  of  Equations  (25),  (30),  and  (31)  must  be  divided  by  | Ay(0,0) | ^ or, 
equivalently,  by  either 

2 2 

/|u(t) | 2 dt  or  /|U(f) | 2 df  (32) 

As  a consequence  of  the  Uncertainty  Principle,  Equation  (19),  the  reverber- 
ation function  Gy(v)  must  also  satisfy  a waveform  - invariant  relation, 

/ Gu(v)  dv  - 1 (33) 
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for  all  normalized  waveforms,  u(t) . From  Equation  (30)  It  is  noted  that  Gu(v)  de- 
pends only  upon  the  spectral  distribution  of  the  transmitted  energy,  |u(f)|^.  To 
get  a feeling  for  the  effect  of  the  signal  bandwidth  on  the  function  Gu(v),  con- 
sider a class  of  signals  whose  spectral  energy  distribution  is  given  by, 


|U(f) |2 


. «f 
stn  1 
nf 

B 


(34) 


This  has  the  same  form  as  the  spectral  energy  distribution  of  a rectangular 
pulse  of  duration  l/B.  It  is  assumed  that  the  waveform  duration  is  T seconds,  and 
the  dependence  upon  T shall  be  indicated  by  writing  Equation  (34)  as 


IV£> 


05) 


where  D - BT  (bandwidth-time  product) . 

When  D ■ 1,  the  waveform  will  be  a simple  pulse,  T seconds  long.  For 
larger  values  of  D,  the  waveform  is  modulated  or  coded  with  0 degrees  of  freedom. 
The  reverberation  functions  corresponding  to  Equation  (35)  are  easily  shown  to  be 


< 


1 ■ 
2a  vT 
0 


T 

) 


sin 


2a  vT 

D 


2a  vT 
D 


(36) 


These  functions  are  plotted  in  Figure  3,  against  the  abscissa,  vT,  which 
is  the  doppler  shift  of  the  target  with  respect  to  the  reverberation,  expressed  in 
units  of  1/T,  the  doppler  resolution  of  the  waveform.  For  low  values  of  doppler 
shift,  the  reverberation  functions  decrease  directly  with  the  dimensionality,  D,  of 
the  waveform.  This  is  to  be  expected  since,  if  the  target  is  not  resolvable  from 
the  reverberation  in  doppler,  the  signal -to-reverberatlon  ratio  can  only  be  improved 
by  increasing  the  range  resolution  and,  therefore,  the  bandwidth  of  the  waveform.  On 
the  other  hand,  for  large  doppler  shifts,  Gp(v)  at  first  Increases  with  increasing  0 
and  then  decreases.  The  bandwidth-time  product,  0,  required  to  equal  the  performance 
of  the  simple  pulse  for  a large  value  of  vT  is  given  approximately  by 


D «* 


(vT)2 


(37) 
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If  this  value  of  BT  can  be  achieved  for  the  largest  doppler  shift  of  in- 
terest, then  it  is  clear  that  a large  BT  waveform  is  the  best  choice  for  combatting 
reverberation.  This  is  fortunate,  since  a large  BT  waveform  is  also  desirable  for 
obtaining  detailed  range  information. 

C . RANDOM  WAVEFORMS 

In  the  preceding  discussions  it  has  been  seen  that  waveforms  having  high 
dimensionality  (bandwidth-time)  products  are  very  desirable  for  target  classifica- 
tion. Wide  bandwidth  alone,  however,  does  not  insure  that  a waveform  will  have  a 
desirable  form  of  ambiguity  function.  The  linear-sweep  FM  waveform  has  a wide  band- 
width and  yet  its  ambiguity  function  indicates  poor  resolution  of  targets  separated 
in  both  range  and  doppler  simultaneously.  One  approach  to  obtaining  a "thumb-tack" 
ambiguity  function  is  to  use  a random  waveform,  that  is,  a waveform  produced  by  a 
random  or  pseudo-random  process. 

In  order  to  discuss  and  describe  the  properties  of  such  waveforms  in  gen- 
eral, statistical  language  must  be  used.  The  exact  ambiguity  function  of  an  unknown 
random  waveform  cannot  be  determined,  but,  if  the  statistics  of  the  process  that  de- 
termined the  waveform  are  known,  the  average  or  expected  value  of  the  ambiguity 
function  can  be  found. 

1.  Gaussian  Noise  Waveforms 

Consider  first  the  class  of  waveforms  that  are  formed  by  gating  the 
output  of  a random  Gaussian  noise  generator  as  shown  in  Figure  4. 

The  noise  at  the  input  to  the  gate  (see  Figure  4)  has  a complex  enve- 
lope, n(t),  whose  power  spectrum,  shaped  by  the  band-pass  filter,  will  be  N0|H(f)|^. 
The  autocorrelation  function  of  n(t)  can  be  found  using  the  Wiener -Khint chin  theorem 
which  gives 

R(x)  - B n(t)  n*(t-x)  (38) 

- N Q / |H(f) | 2 eJ  2n  fT  df 

It  is  assumed  that  the  amplitude  of  H(f)  is  adjusted  so  that  R(o)  “ 1. 
This  requires  that 

No  / |H(f)|2  df  - 1 (39) 
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The  waveform  at  the  gate  output  is  the  transmitted  waveform,  u(t). 


which  may  be  described  as 


> for  I I < T/2 


u(t)  - 


0 , for  |t|  > T/2 

The  ambiguity  function  of  u(t)  is  defined  to  be 
T/2-t 


Au(t,v)  -f  u(t)  u*(t  + t)  e'^  2n  vt 


dt 


-T/2 


and  therefore,  the  expectation  or  statistical  average  of  A(t,v)  will  be 

T/2-t 


E- 


j^Au(T,v)|-  •J'  E-^u(t)  u*(t  + t)J-  e”^  2n  vt  dt 

-T/2 

Using  Equations  (38)  and  (40) 

E*(t),  for  | t | and  |t  + t|  < T/2 


E-<  u(t)  u*(t  + t) 


0 , otherwise 

which  enables  evaluation  of  Equation  (42) ; the  result  is 

T - | t | sin  fl  v (T  - | t | ) 


E 


-|^Au(T,v)j-  - R*(t)  n VT  T 

- R*(t)  At(t,v) 


n v (T  - | x | ) 


(40) 


(41) 


(42) 


(43) 


(44) 


where  A^(t,v)  is  the  ambiguity  function  of  a simple  pulse,  T seconds  long. 

Assume  now  that  the  bandwidth  of  the  filter  H(f)  is  much  broader  than 
1/T,  so  that  u(t)  has  a high  bandwidth-time  product.  In  this  case,  the  width  of  the 
correlation  function  R(t),  will  be  very  small  compared  to  T.  As  a result,  the  range 
dependence  and  range  resolution  of  the  expected  ambiguity  function,  E Ay  , will  be 
determined  essentially  by  the  shape  and  width  of  R(t) . The  doppler  dependence  of 
E j^AyJ  is  determined  entirely  by  the  doppler  dependence  of  A^(t,v),  which  has  a 
doppler  resolution  of  1/T. 
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Although  the  expected  or  average  ambiguity  function  calculated  above 
gives  some  Indications  of  the  resolution  performance  of  a randomly  selected  wave- 
form, It  is  an  incomplete  description  and  in  some  respects  misleading.  It  is  known 
that  the  ambiguity  function  of  any  particular  waveform,  u(t),  must  satisfy  the 
Uncertainty  Principle.  The  average  ambiguity  function  given  by  Equation  (44),  how- 
ever, does  not  satisfy  the  Uncertainty  Principle,  since 


II  |R(t)|  |At(t,v)|2  dx  dv  < II  |At(x,v)|2  dx  dv  - 1 


(45) 


A statistical  measure  of  the  resolution  performance  of  a randomly  se- 
lected waveform,  u(t),  that  does  satisfy  the  Uncertainty  Principle  can  be  obtained 
by  finding  the  expected  value  or  statistical  average  of  ^(Xjv))2.  This  makes  more 
sense  physically  since  the  signals  that  are  usually  of  concern  are  those  that  add 
power-wise  rather  than  voltage-wise. 

The  expectation  of  |Ay|2  can  be  derived  as  follows 
T/2-x 

f . ■ i.  . . 1 

dt  ds  (46) 

_ _ i i 

-T/2 


e||Au(x,v)|2|  mJJ  s|u(t)  U*(t  + |x|)  u*(s)  u(s  + |x|)}  e‘j2nv(t's) 


T/2- | x | 

^ 1 1 'R<T>^2  + IR(t‘8)|2  e'j2nv(t's)  dt  ds 

izhl 

2 

|R(t)|2  Ia^XjV)!2  + ~~  ff  | R(t-s)  | 2 d-J2nv(t-s)  dt  ds 

T-lxl 


(47) 


(48) 


T-  x 


- | R(x)  | 2 |AT(x,v)|2  | R(x)  |2  (T- 1 x | - |x| ) e‘J2nvx  dx 


- (T-  x ) 


(49) 


For  waveforms  having  large  BT  products,  |R(x)|^  will  have  a small 
spread  in  "x"  compared  to  T-|x|  (for  most  values  of  x),  and  hence  the  following  ap- 
proximation to  Equation  (49)  may  be  obtained  for  use  of  Equations  (38)  and  (39) . 


E{lV^v)|2j*.  I*<t)|2  IVx,'')!2  + "V  - r(v) 


(50) 
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where 


T(v)  » 


/ | H(f ) | 2 | H(f+v) | 2 df 
/ | H(f) | 2 df  2 


Assuming  further  that  | H(f ) | 2 is  approximately  rectangular  in  shape 
with  a bandwidth  of  B cps,  Equation  (50)  can  be  simplified  to 

E ||A(t,v)|^«  |R(T)|2  ^(T^)!2  + ■—  1 - 1 - -T-  (52) 

The  first  term  in  Equation  (52)  is  the  square  of  the  value  obtained 
previously  for  the  average  value  of  Ay.  It  represents  the  central  spike  in  the 
ambiguity  function.  Its  range  resolution  corresponds  to  the  bandwidth  of  u(t)  and 
its  doppler  resolution  to  the  dime  duration  of  u(t).  The  second  term  is  the  unde- 
sirable part  of  the  ambiguity  function,  that  is,  the  side  lobes  or  self-noise  of 
the  waveform.  In  the  t,v  plane  it  extends  over  ±T  in  t and  ±B  in  v with  a peak 
height  of  1/BT  and  a total  volume  of  unity. 

The  performance  of  the  randomly  selected  waveform  against  reverbera- 
tion is  considered  next.  It  has  been  seen  that  for  an  explicitly  defined  waveform 
this  is  determined  by  the  reverberation  function 


G (v) 
u 


* 

-/  |Au(t,v)|' 


For  a randomly  selected  waveform,  the  average  performance  in  reverber- 
ation is  determined  by  the  expected  value  of  Gy 


i 

£„<*>}  •/  * I } ■*' 


Using  Equation  (49),  this  becomes 


& 

^u(v>!  mJ  1*00  1 2 Vt,v)|2  dT 


-V  f f |»(*)|2  (T-|t|-1x|)  e"j2nvx  dx  dt 

T -T  -(T-|t|)  (J 


I 
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1 


M -j C 


, sin  nv(T-lxl) 

i R(x)  \ 5 dT 

(nvT)Z 


+ -y* y*  | R(t)  1 2 (T- 1 t I ) 2 e'*2™  dx 
T -T 


(56) 


Assuming^  as  before,  that  u(t)  is  a large  BT  product  waveform  so  that 
| R(t) | 2 is  a delta-like  function  compared  to  the  other  terms  in  Equation  (56),  the 
following  approximation  is  obtained 


E 


2 

{Gu(v)}  “ StnvTVT'  / lR<T)|2  dT  + / iR<T>|2  e‘j2nvT  dx  (57) 


Both  integrals  in  Equation  (57)  can  be  expressed  in  terms  of  the  func- 
tion r(v)  defined  by  Equation  (51),  so  that 


E 


{v*)} 


8 in  itvT 
nvT 


r(o)  + r(v) 


(58) 


If  the  power  spectrum  of  the  noise  process  is  approximately  rectangular 
in  shape  with  a bandwidth  of  B cps,  the  expected  value  of  Gy  becomes 

,2 


E 


{v->} 


T 

BT 


J\  . _l£_  + r_§ln,.^T_-| Z ' 

\L  3T  + L nvT  1 J 


(59) 


It  is  interesting  to  compare  this  result  with  the  reverberation  func- 
tion of  a "definite"  waveform  that  is  T seconds  long  and  has  a rectangular  spectral 
distribution  of  energy  of  the  form 


|U(f)|‘ 


’] 

■ 


(60) 


"l/B  for  |f|  < B/2 
0 for  |f|  > B/2 
Substituting  Equation  (60)  in  Equation  (27) 

°u(v)  ■ -$r  [l  • i*-]  «■*> 

Comparing  this  with  Equation  (59),  it  is  seen  that  the  random  waveform 
has  approximately  3 db  poorer  performance  against  reverberation  for  values  of  vT  < j; 
for  larger  values  of  vT  the  random  waveform  does  as  well  as  the  "definite"  waveform. 
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. Random  Binary-Phase  Modulation 

The  random  waveforms  considered  in  the  foregoing  discussion  were  pro- 
duced by  gating  T seconds  of  band-pass  Gaussian  noise.  As  a result,  these  waveforms 
were  modulated  in  both  amplitude  and  phase  within  the  T-second  gate  period.  It  is 
very  desirable,  however,  to  use  waveforms  that  have  a constant  amplitude  during  the 
gate  period  and  are  coded  or  modulated  in  phase  only.  Transmitter  efficiency  is 
usually  higher  with  constant  amplitude  waveforms  and,  more  important,  maximum  energy 
is  obtained  for  a given  pulse  length  since  cavitation  limits  the  peak  power  delivered 
to  the  water. 

A simple  example  of  a randomly,  phase -modulated  waveform  is  a carrier 
whose  phase  is  changed  by  either  0°  or  180°  at  periodic  intervals.  Such  a waveform 
could  be  generated  by  multiplying  a sinusoid  carrier  by  a random  square  wave,  such 
as  .shown  in  Figure  5.  If  the  waveform  of  Figure  5 la  called  u(t),  the  complex  rep- 
resentation of  the  complete  signal  would  be  u(t)  exp  J2xfct.  To  get  a mathematically 
tractable  formulation  of  u(t),  note  that  u(t)  may  be  considered  to  be  the  sum  of  a 
sequence  of  identical  pulses,  each  of  duration  X,  with  randomly  chosen  algebraic 
signs,  that  Is 

N-l 

u(t)  - ^ 7^  p (t  - kX)  (62) 

k=0 


It  shall  be  assumed  that  the  random  variables,  xk,  are  independent, 


so  that 


r ^ 0 for  k i l 

7 " £or  k .i 
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As  before  Che  expected  values  of  the  ambiguity  function,  A^TjV)  shall 
be  calculated,  and  the  expected  value  of  the  reverberation  function,  Gu(v) . 

By  definition 

Au(t,v)  ■ / u(t)  u*(t  + t)  e ^ 2n  vt  dt  (66) 

Substituting  from  Equation  (62) 


N-l 


Au(t,v)  *y  y *k  x^  / p(t  * kX)  p(t  + T -ix)  e 2n  vt  dt  (67) 

k,  f=6 

This  can  be  written  very  conveniently  in  terms  of  the  ambiguity  func- 
tion of  the  basic  pulse  p(t) . Thus,  letting 

-1  2*  vt  in  vX  8ln  «v(X-M) 

+ t)  e J VC  dt  - eJ  vA  (68) 

nv(X- | x | ) 


Ap(t,v)  - f p(t)  p(t  + t) 


Equation  (67)  can  be  written  as 

N-l 


x,  x,  e'j  2n  k vX  A 


j(k-jt)  X + T,vj 


(69) 


The  expected  value  of  AyOr^v)  is  easily  found  using  Equation  (65) 

N-l 

£ E{*k  2*  k VX  X + T,vJ  (70) 

• -f  I •'*  2”  k VX  Ap(l'v> 

k«0 


- e-J(N-D«vX  Jte  lgy?,.,  k ( j 
N sin  nvX  p ' 

The  expected  value  of  ^1^  is  somewhat  more  difficult  to  obtain.  Starting  with 

k,i,r,s-0 


x A 


[(k.|)X.7,v]V[(r..)X  + „v] 
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and  noting  that  E -j  x xr  xg  >•  is  either  zero  or  1/N2  It  will  equal  1/N2  if,  and 
only  if  (a)  kmJL  ana  r«*s,  or  (b)  k«r  and  l»s,  or  (c)  k»s  and  lm r.  These  three  con- 
ditions are  not  mutually  exclusive  since  they  contain  the  common  condition 
(d)  k«i»r»s . 

can  be  obtained  by  adding  the  contributions  to  the  quadruple 
sum  on  the  right-hand  side  of  Equation  (71)  under  conditions  (a),  (b) , and  (c),  and 
then  subtracting  twice  the  contribution  to  the  quadruple  sum  under  condition  (d) . 

The  contributions  for  each  of  the  four  cases  are 

la)  7 1 I e-J2”vX<k-r>  i*p(t,v»i2  - i*p(,,v)i2  [s1;^]2  <72> 


k,  r-0 


N-l  N-l 

<b)  i 1 1 i*p  [lk-l>  x + T’vJi2  ■ I 


k,  £.-0 


k— (N-l) 


~2~  |Ap  kX+x,v  | (73) 


^ Z Z e’J2flvx<k’t)  {AP[(k'4)  x+  T'v]}iv  (l~k)  x + t2v]]><74) 


k,  i-0 


t X iAP(T'v)i2  mjr  iAP<T'v>i2 

N I A 


The  summation  for  case  (c)  reduces  to  1/N  | Ap (t , v) | ^ since 
Ap(x, v)  ■ 0 for  | t | > X 

Assembling  the  above  results 


*{[VT.»i2}-[-£1:1”n!£x]  IV’^I2  - i iy^>l: 

T N'|k|  , ,2 

+ 2^  ~^2 lAp(^x  + T,v)r 


k-(N-l) 


The  first  term  on  the  ’■ight-hand  side  represents  the  "central  spike"  of 
the  ambiguity  function  Its  width  in  doppler,  determined  by  the  first  factor,  is 
1/T  where  T » NX,  is  the  total  duration  of  the  waveform.  The  width  in  the  range 
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dimension  is  determined  by  the  second  factor  | Ap | ^ and,  therefore,  is  approximately 
X or  T/N.  Note  that  N is  a measure  of  the  dimensionality  or  bandwidth-time  product 
of  the  waveform.  The  remaining  two  terms  in  Equation  (77)  represent  the  ambiguity 
side  lobes  or  self-noise  of  the  waveform.  The  side  lobes  have  a maximum  value  of 
approximately  1/N  at  t,v  ■ ± X,0  and  decrease  linearly  with  increasing  |t|. 

The  expected  reverberation  function,  E-^G^,  18  now  evaluated.  Apply- 
ing the  expectation  operator  to  Equation  (25) 

e{cu(v)J  - / E||Au(T,v)|2JdT  (78) 

and  then  using  the  expression  for  the  Integrand  given  in  Equation  (77),  it  is  found 
that 


«{c,(v>}-  Gp<v){-Sjl-  + (79) 

where 

yv>-nAp<T,v>i2dT.-^-{1-  *>■#»}  (eo, 

Substituting  this  expression  in  Equation  (79),  and  setting  T - NX 

* {«.<*>}-¥  pfej 

Except  for  the  last  factor,  this  result  is  identical  to  the  result  ob- 
tained in  Equation  (36)  for  the  reverberation  function  of  a "definite"  waveform 
having  the  energy  spectral  distribution  of  a rectangular  pulse  and  a dimensionality 
of  N.  The  second  factor,  however,  Indicates  that  the  average  reverberation  function 
for  the  random  waveform  is  approximately  3 db  larger  for  values  of  vT  < 1/2.  This 
is  similar  to  the  result  obtained  for  Gaussian  random  waveforms. 


i 2«vT 
sin  - .1  •• 

„ . 2 

. N 

N-l  . isln  NkvX  j 

flail 

N In  sin  wvX 

L N J J 

(81) 
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